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Late evolution of the X-ray afterglow of GRB 030329 
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Abstract. The X-ray afterglow of the Gamma-Ray Burst GRB 030329 , associated to SN2003dh at ^=0.1685, has 
been observed with XMM-Newton 258 days after the burst explosion. A source with flux of (6.2±2.3) x 10"^'' erg 
cm~^ s~^ (0.5-2 keV) has been detected at the GRB position. This measurement, together with a re-analysis of 
the previous X-ray observations, indicates a flattening of the X-ray light curve ~40 days after the burst. This is in 
remarkable agreement with the scenario invoking the presence of two jets with different opening angles. The wider 
jet should be responsible for the observed flattening due to its transition into the non-relativistic Sedov-Taylor 
phase. 
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1. Introduction 

Multifrequency studies of Gamma-Ray Burst (GRB) af- 
terglows at late times ( ^ 10 days) can be used to test 
some predictions of the standard, as well as alterna- 
tive, models. This has been done so far for a handful of 
bursts (e.g., Frail, Waxman & Kulkarni 2000; Frail et al. 
2004). Thanks to its high flux and low redshift (z=0.1685, 
Greiner et al. 2003; Caldwell et al. 2003), GRB 030329 is 
an ideal target for such kind of investigations. 

In fact GRB 030329 had a very large fluence of ^lO'^ 
erg cm-2 (30-400 keV, Ricker et al. 2003), in the top 1% 
of all detected GRBs, and its optical transient had mag- 
nitude 13 one hour after the explosion (Peterson & Price 
2003; Torii 2003). GRB 030329 is also the first GRB unam- 
biguously associated with a supernova (SN2003dh, Stanek 
et al. 2003; Hjorth et al. 2003). Systematic multiwave- 
length monitoring of its afterglow (e.g., Lipkin et al. 2004; 
Matheson et al. 2003; Tiengo et al. 2003; Sheth et al. 2003) 
is yielding important information for the understanding of 
the jet structure, GRB energetics and circumburst envi- 
ronment. 

Here we report on a new XMM-Newton observation 
performed 258 days after the burst in order to study the 
late evolution of the X-ray afterglow. We also re-analysed 
the previous observations to derive in a consistent way all 
the X-ray fluxes in the 0.5-2 keV energy range. This is the 
range where the effective area of the XMM-Newton EPIC 
instrument is highest and the corresponding fluxes are less 
affected by uncertainties in the afterglow spectral shape. 



2. Data analysis and results 

XMM-Newton observed the position of GRB 030329 start- 
ing on December 12, 2003 at 7:30 UT, for an observation 
length of ~24 hours. After excluding the time intervals af- 
fected by high particle background, the net exposure times 
were 68 and 73 ks, respectively in the PN and MOS cam- 
eras of the EPIC instrument (Striider et al. 2001; Turner 
et al. 2001). All the cameras operated in Full Frame mode 
and with the thin optical blocking filter. The data were 
processed using SAS version 5.4.1. 

A faint source at the GRB position is visible, albeit 
barely, in the 0.5-10 keV images. The previous X-ray ob- 
servations (Tiengo et al. 2003) indicated a rather soft spec- 
trum and a low absorption, which, combined with the high 
effective area of the PN camera at very low energies, re- 
sults in a significant fraction of the source counts in the 
0.2-0.5 keV band. To increase the signal to noise ratio 
we therefore analyzed the data over the larger 0.2-10 keV 
range. This required an additional cleaning procedure^ to 
reduce the instrumental noise, which in the PN detector 
increases very rapidly below '^0.4 keV. The resulting im- 
age is shown in Fig.l, where the small solid circle indicates 
the GRB 030329 position. 

The estimate of the afterglow flux requires particular 
care, due to the proximity (~30") of an X-ray emitting 
AGN, which in this observation was much brighter than 



^ we processed the PN data using an experimental version 
of the task "epchain" included in a development version of the 
SAS software 
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the GRB afterglow. For this reason we followed different 
methods to estimate the background and the fraction of 
AGN counts falling in the source extraction region. For 
the afterglow counts we used a circular region with radius 
of 15" centered at the accurately known GRB position 
(Taylor et al. 2003). The total number of counts (i.e. back- 
ground, AGN contamination, and afterglow) contained in 
this region is 148 for the PN and 174 for the sum of the 
two MOS. 

With the first method, we estimated separately the 
background and the AGN contamination. The former was 
derived from a source free region (a circle of radius 40") in 
the same CCD chip, while our estimate of the latter was 
based on the well known model of the instrumental Point 
Spread Function (PSF, Ghizzardi 2002). We measured 
the net AGN counts within 15" from the AGN position 
(572±26 in the PN and 456±25 in the two MOS summed). 
The PSF model predicts that a number of counts cor- 
responding to 5% and 4% of these values (for the PN 
and MOS cameras, respectively) will fall in the 15" ra- 
dius GRB extraction region. Thus we finally obtain an 
estimate of 50±19 net afterglow counts in the 0.2-10 keV 
energy band (33±13 in the PN camera alone). 

An alternative approach consists in directly measuring 
the background and AGN contamination based on an ap- 
propriate extraction region close to the afterglow position, 
as was done in Tiengo et al. (2003) In this way we obtain 



54±19 net afterglow counts in the 15" GRB extraction 
region (see Fig. 1). 

As a further test we developed a maximum likelihood 
analysis program and used it to model the region of the 
0.2-10 keV PN image containing GRB 030329 and the 
AGN. We used the appropriate PSF and let the positions 
and intensities of the two sources, as well as the back- 
ground level, as free parameters. This yielded an intensity 
normalization for the afterglow corresponding to 34 PN 
counts within a radius of 15". The source positions, the 
background level and the AGN intensity derived in this 
way were also consistent with the values found with the 
methods described above. 

In conclusion, since all methods gave consistent results, 
in the following we adopt the count estimate obtained with 
the first method, which corresponds to a l.la detection. 
Since the small statistics does not permit to carry out a 
detailed spectral analysis, a spectral model must be as- 
sumed to convert count rates to physical flux units. The 
spectra measured in the previous X-ray observations were 
well fitted by a power law. In the first Rossi-XTE ob- 
servation the photon index could be tightly constrained 
at r=2.17lQQ3. '^^'^ following observations, performed by 
XMM-Newton , gave a strong upper limit on the neu- 
tral hydrogen column density, which cannot differ sub- 
stantially from the Galactic value of Nh=2x10^*' cm~^. 
We have then verified that an absorbed power law with 
r=2.17 and = 2 x 10^° cm~^ gives an acceptable fit to 
all the X-ray spectra of the GRB 030329 afterglow mea- 
sured up to now. We thus assume such a spectrum, which 
yields an observed flux (i.e. not corrected for the absorp- 



Time since burst 


Satellite 


0.5-2 keV flux 


(days) 




(10^^^ erg cm"^ s"^) 


0.2065 


Rossi-XTE 


153±15 


0.2125 


Rossi-XTE 


146±15 


0.2185 


Rossi-XTE 


162±15 


0.2515 


Rossi-XTE 


134±18 


1.2762 


Rossi-XTE 


9.3±3 


37.3 


XMM-Newton 


0.0143±0.0007 


60.9 


XMM-Newton 


0.0079±0.0005 


258.3 


XMM-Newton 


0.00062±0.00023 



Table 1. Fluxes of the afterglow of GRB 030329 in the 
0.5-2 keV energy range. 



tion) of (6.2±2.3)x 10^16 erg cm^^ g- 
range. 



in the 0.5-2 keV 



3. The X-ray afterglow long term evolution 

To derive the long term flux history of the afterglow in 
a consistent way, we have re-analyzed the two XMM- 
Newton observations of May 2003 using the procedures 
described above. In particular, we reprocessed the PN data 
with the new cleaning algorithm and verified that the dif- 
ferent ways of estimating the AGN contamination gave 
consistent results. 

As for the December observation, we finally adopt the 
values based on the first method: 501±25 and 400±25 
net afterglow counts for the first and second observa- 
tion, respectively. The light curve from the three XMM- 
Newton observations is well fit by a power law decay with 
index (5=1.45±0.2 (Fig. 2). 

For a power law spectrum with r=2.17 and A'h — 
2 X 10^*^ cm~^, the above count rates correspond to 0.5-2 
keV observed fluxes of (143±7) x 10"^^ and (79±5) x lO'^^ 
erg cm~^ s~^, respectively ^. These flux estimates do not 
change significantly if the spectral parameters are allowed 
to vary within their 90% confidence intervals (Tiengo et 
al. 2003). 

To draw a complete X-ray light curve containing both 
the XMM-Newton and Rossi-XTE data, an extrapolation 
of the X-ray spectrum is unavoidable. In fact, most of the 
counts detected by EPIC have energies below 2 keV, while 
the PGA instrument of the Rossi-XTE satellite is not sen- 
sitive to such low energy photons. We therefore extracted 
the 0.5-2 keV fluxes of the Rossi-XTE observations from 
the best fit absorbed power law model to each PCA spec- 
trum, obtaining the values reported in Table 1. The cor- 
responding errors take into account the uncertainties on 
the spectral shape. The resulting light curve is shown in 
Figure 3. Fitting it with a single power law yields a time 
decay with (5=1.78±0.02 and reduced of 2.6. 



^ the corresponding values in the 0.2-10 
((351±17)xl0"^'' and (194±12) x 10"^'' 



keV 

-2 



range 

erg cm " s ^) are 
slightly different, but consistent within the uncertainties, from 
those reported in |Tiengo et al. (2003)| which were derived from 
spectral fitting 
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Fig. 1. PN image (0.2-10 keV) of the field of GRB 030329 
on December 12-13, 2003. The extraction region for the 
afterglow (solid circle) and one of the regions used for 
the background subtraction are shown. The background 
region is an annulus centered at the position of the nearby 
AGN with inner and outer radii of 15" and 45", with the 
exclusion of the 15" circle around the GRB and other two 
weak sources. 



4. Discussion 

As shown above, a single power law decay gives a 
marginally acceptable (1.6% c.l.) description of the X-ray 
afterglow from 5 hours to almost nine months from the 
burst. 



The data presented in Tiengo et al. (2003) 



were con- 
sistent with a break in the X-ray light curve at t~0.45 
days. Similarly to the optical data available at that time, 
the decay slope changed from ~0.9 to ~1.9. Note that in 
Tiengo et al. (2003) we used fluxes in the 2-10 keV en- 
ergy band, which was optimal to represent the observed 
flux in the PCA range, but implied a substantial extrap- 
olation of the poorly constrained XMM-Newton spectra. 
This resulted in large errors on the XMM-Newton 2-10 
keV fluxes. Due to the different energy range used here 
(0.5-2 keV), the points at t=37 and 61 days have smaller 
errors. They are not consistent with a slope of 1.9 but are 
instead fitted by a fiatter power law (see Fig. 2). 

In the following we give a more complex interpretation 
of the X-ray data drive n by the two-jet model proposed 
by Berger et al. (2003) Although this is not the unique 



way to fit the X-ray light curve, it has the advantage of 
explaining the complete X-ray data set presented here as 
a natural consequence of an already published interpreta- 
tion of the light curves at optical and radio wavelengths. 

On the basis of multifrequency radio data, 
Berger et al. (2003)| suggested that GRB 030329 is 
characterized by a two-component jet, with different 
opening angles (5 and 17 degrees). The wide jet has 
an intrinsic kinetic energy, corrected for collimation, of 
E = 2.5 X 10^" erg, while the narrow jet emits an intrinsic 



50 100 200 

Time since Burst (days) 

Fig. 2. Count rate of the GRB 030329 afterglow in the 
0.2-10 keV energy band as measured by the EPIC in- 
strument (the plotted count rate is from the sum of the 
PN, MOSl and M0S2 net counts within an extraction 
region of 15"). The line is the best-fit power law model 
{S = 1.45±0.2). 



7-ray energy of i?^ = 5 x 10^^ erg. Both jets are assumed 
to move in an interstellar medium with a constant density 
of n = 1.8 cm^'^. The initial part of the optical and X-ray 
afterglow should originate from the narrow jet, while the 
wide jet should be responsible for the radio flux, the late 
X-ray afterglow, and it also explains the resurgence of 
the optical flux after ~1 day (Berger et al. 2003; Sheth et 
al. 2003; Lipkin et al. 2004). 

The isotropic equivalent energies of the two jets are 
6.6 X 10^^ erg (narrow jet, assuming a kinetic to gamma- 
rays conversion factor of ?7^=0.2) and 5.7 x 10*^^ erg (wide 
jet). As a consequence, the wide jet should become non- 
relativistic earlier than the narrow jet. The transition to 
the non-relativistic Sedov-Taylor phase manifests itself 
through a change in the decaying slope of the flux, which 
becomes F{u,t) oc where Snr = (15p - 21)/10 in 

the case of constant density medium, and p is the slope of 
the injected relativistic electron distribution. In our case, 
for any reasonable value of p, the transition should give a 
flattening of the light curve, if it occurs after the jet break. 

The time tjyii of such a transition is not unambiguously 
predicted theoretically, but for GRB 970508 (z=0.835) it 
was seen to occur at t^j^ ^ 100 days (Frail, Waxman & 
Kulkarni 2000). Since in that case an isotropic equivalent 
total kinetic energy of E' = 5 x 10^^ erg and n = 0.5 cm~'^ 
were estimated, we can assume 



tNR ^ 55(1 



1/3 



days. 



(1) 



This value is also consistent with theoretical expectations 
(see e.g. Dai & Lu 1999; Livio & Waxman 2000; Huang, 
Dai & Lu 1998) and also with the interpretation given 
by [Berger (2003)] of the radio light curve of GRB 980703. 
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With the parameters derived by 'Berger et al. (2003) for 
the wide jet component of GRB 030329, we obtain 



tNR. ~ 44 



E 



5.7 X lO^ierg 



1/3 



1.8 

n 



1/3 



days, 



(2) 



which is close to the time of the first XMM~Newton ob- 
servation. The flat time decay {S = 1.45 ± 0.2) can thus 
be explained by the wide jet component which becomes 
non-relativistic. This value of Snr gives a slope of the 
electron distribution p = 2.37 ± 0.13, corresponding to an 
X-ray spectrum with T = {p/2 + 1) = 2.18 ± 0.07. The 
XMM-Newton spectra are consistent with this value^. 

At these late times, the cooling frequency is already 
below the optical band. Therefore, the optical and the 
X-ray light curves should have the same decay slope. 
After accounting for the contribution from the super- 
nova, Berger et al. (2003)| and'Lipkin et al. (2004) derived 
Sopt ~ 2.35 between day 10 and day 40 (i.e. after the jet 
break of the wide jet). This is in agreement with the ex- 
pected value Sopt = P l |Sari, Piran fc Halpern T999| ). Note 
that the optical light curve after 40 days is dominated by 
the supernova light, which masks the predicted flattening. 

Finally, we would like to caution the reader that there 
may exist other possible interpretations of the observed 
flattening of the X-ray light-curve. In the cannonball sce- 
nario (e.g. Dado, Dar & De Rujula 2002) a flattening of 
the light curve can be associated to the afterglow contri- 
bution of a new cannonball: indeed, in Dado, Dar & De 
Rujula (2004) the optical afterglow of GRB 030329 is ex- 
plained with two cannonballs (plus the SN contribution) 
up to 70 days after the trigger. A flattening of the X-ray 
light-curve could then be possible if a third cannonball 
starts to contibute in the X-ray band (but not in the op- 
tical, which is dominated by the SN light). However, the 
number of cannonballs should be equal to the number of 
main pulses in the prompt 7-ray emission, which are only 
two in GRB 030329. 

In the dyadosphere model (e.g. Ruffini et al. 2003), a 
flattening of the light-curve is predicted, as in the more 
conventional flreball scenario, when a transition to the 
Newtonian expansion occurs, but other features specific 
to GRB 030329 (such as the break at ~ 0.5 days) would 
remain unexplained. 

5. Conclusions 

The combination of the exceptional brightness of GRB 
030329 and the high sensitivity of XMM-Newton has al- 
lowed us to study a GRB X-ray afterglow up to very late 
time, leading to a ;S 3(t detection of the afterglow 258 
days after the burst. 

Although other interpretations are possible, the 
XMM-Newton data at t =37, 61 and 258 days are con- 
sistent with the X-ray light curve being dominated by 



^ the coincidence with the spectral slope derived by the first 
7?ossi-XTiJ observation is, in this scenario, a pure coincidence, 
since the early X-ray emission is due to the narrow jet. 
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Fig. 3. Optical (Lipkin et al. 2004) and X-ray [0.5-2 keV] 
light curves of the GRB 030329 afterglow. The two solid 
lines (one for the X-ray and one for the optical) corre- 
spond to the prediction of the narrow-|-wide jet scenario, 
modified at late times to account for the non-relativistic 
phase of the wide jet. Dotted lines indicate the contribu- 
tion of the wide jet only, the short-dashed line corresponds 
to the optical contribution of the supernova, and the long- 
dashed line corresponds to the optical contribution of the 
narrow jet, with a break at 12 hours. The jet-break time 
for the wide component is assumed to occur 10 days after 
the burst. 



the wide jet component envisaged by Berger et al. (2003), 
Lipkin et al. (2004) and Sheth et al. (2003). In particular, 
and more importantly, we have evidence that this compo- 
nent became non-relativistic at ~40 days after the GRB 
explosion, in remarkable agreement with what predicted 
assuming the parameters (kinetic energy and density of 
the external medium) derived by the above authors by 
fitting the radio and optical light curves. The immedi- 
ate prediction of this interpretation is that also the ra- 
dio light curve should show a flattening from 6 — 2 to 
Snr^ ^ 1-5 after t ~40 days. The publicly available ra- 
dio data at these late times have large uncertainties and 
cannot exclude such a behavior. Only further analysis of 
more radio data could settle this issue. The lack of a flat- 
tening in the radio light curve would imply that the fire- 
ball has not yet reached the non-relativistic phase, due 
to a larger total kinetic energy of the wide jet and/or a 
smaller density for the interstellar medium. This would 
then require a different reason for the observed flattening 
in the X-ray light curve, e.g. a larger supernova contri- 
bution in the X-ray band. This would imply an X-ray 
luminosity at least one order of magnitude higher than 



A. Tiengo et al.: Late evolution of the X-ray afterglow of GRB 030329 



5 



that observed in other type Ic and II supernovae (see Fig. 
3 of IKouvehotou et al. 2004|l . 
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